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Abstract. We use the quasiclassical theory of superconductivity to calculate the re- 
sponse of the currents in the core of a vortex to an alternating electric field. We consider 
relatively clean superconductors with a mean free path, £, larger than the coherence 
length, £o- The current response of the core is dominated by the bound states of Caroli, 
de Gennes and Matricon, and differs significantly from the response obtained from the 
model of a "normal core" introduced by Bardeen and Stephen. The model of Bardeen 
and Stephen describes the limit of a dirty core (£ <C £o), and fails for a clean supercon- 
ductor. The response of the bound states includes non-dissipative acceleration of the 
charge carriers as well as dissipative currents. Dissipation dominates at low frequencies 
(hui ^ -5A), and lZetj(u>) greatly exceeds the normal state conductivity. These dissi- 
pative currents have a dipolar structure superimposed on an oscillating supercurrent 
in the direction of the electric field. We attribute the anomalous dissipation to the 
zero-energy bound state broadened by lifetime effects. 1 

1 Introduction 

An electric field has two principal effects on a superconductor. It changes the 
supercurrents by accelerating the superfluid condensate, and it generates dissipa- 
tion by exciting "normal" quasiparticles. This two- fluid picture of a condensate 
and normal excitations is clearly reflected in the optical conductivity of standard 
s-wave superconductors in the Meissner state (Mattis and Bardeen 1958). The 
conductivity has a superfluid part, 



and a dissipative normal part. At T = dissipation starts at the threshold 
frequency for creating quasiparticles, u> — 2A. At finite temperature the coupling 
to thermally excited quasiparticles leads to dissipation at all frequencies. The 
two-fluid picture no longer holds for type-II superconductors in the vortex state. 
The response to an electric field consists of contributions from gapless 2 vortex 
core excitations (Caroli et al. 1964), and contributions from excitations outside 
the cores, where the quasiparticle spectrum shows the bulk gap. For high-K 
superconductors the electric response at low frequencies of the region outside 
the core can be described very well by the London equations (F. London 1950), 

1 This work was done in collaboration with Dierk Rainer and James A. Sauls. 

2 We ignore the "mini-gap" of size A 2 /Ef predicted by Caroli et al. 1964. 
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i.e. by two-fluid electrodynamics. The response of the core is more complex. 
In the traditional model of a "normal core" (Bardeen and Stephen, 1965) the 
conductivity in the core is that of the normal state, a core = a n . The Bardeen- 
Stephcn model is a plausible approximation for dirty superconductors with a 
mean free path (£) much shorter than the size of the core [l <C £o)- In this 
limit the vortex core excitations of Caroli et al. 1964 may be considered as a 
continuum of normal excitations. 

In 1969 Bardeen et al. published a detailed discussion of the bound states in 
the core of a vortex, and argued that the bound states contribute significantly to 
the circulating supercurrents in the vortex core. This effect was confirmed by re- 
cent self-consistent calculations of free and pinned vortices in clean and medium 
dirty superconductors (Rainer et al. 1996). The authors show that all currents 
in the core (circulating currents as well as superimposed transport currents) are 
carried predominantly by the bound states. This means that the model of a 
normal core needs to be modified for relatively clean (£ ^ £ ) superconductors. 
Furthermore the response of the bound states is expected to show dissipative as 
well as superfluid features. 

In this article we present calculations of the linear response of the current 
in the core of an isolated vortex to an a.c. electric field with frequency below 
the bulk threshold frequency 2 A In section 2 we describe the model as well as 
the basic theoretical concepts and equations of the quasiclassical linear response 
theory of inhomogeneous superconductors. We discuss the importance of vertex 
corrections, charge conservation, and the local charge-neutrality condition. Self- 
consistent calculations of the structure and excitation spectrum of a vortex in 
equilibrium arc presented in section 3. The analysis and discussion emphasizes of 
the effects of impurities on the bound states in the core. Our numerical methods 
are based on a formulation of the quasiclassical response theory in terms of 
"distribution functions" , which are introduced in section 4. In sections 5 and 6 
we present and discuss our results for the current response to a small a.c. electric 
field. 

2 Quasiclassical theory of the electromagnetic response 

Our calculation of the response of a vortex to an a.c. electric field is based on 
the quasiclassical theory of Fermi liquid superconductivity. This theory describes 
phenomena on length scales large compared to the microscopic scales (Bohr ra- 
dius, lattice constant, kj , Thomas-Fermi screening length, etc.) and frequencies 
small compared to the microscopic scales (Fermi energy, plasma frequency, con- 
duction band width, etc.). 3 The a.c. frequencies we consider are of the order of 
the superconducting energy gap, | A |~ T c , or smaller, and the length scales of 
interest are the coherence length, £ = Vf/2nT c , and the penetration depth, A. 
Hence, our theory requires the conditions kf£o 3> 1, and T c /Ef <C 1, where the 

3 We set in this article % = Ub = 1. The charge of an electron is e < . 
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Fermi wavelength, kj 1 , and the Fermi energy, Ef, stand for typical microscopic 
length and energy scales. 

A convenient formulation of the quasiclassical theory for our purposes is in 
terms of the quasiclassical Nambu-Keldysh propagator g(pf, R; e, t), which is a 
4 x 4-matrix in Nambu-Keldysh space, and a function of position R, time t, 
energy e, and Fermi momentum pf. We consider a superconductor with random 
atomic size impurities in a static magnetic held, and an externally applied a.c. 
electric field, E = —-d t 6A. In a compact notation the transport equation for 



this system and the normalization condition read 4 

[(e + -w s • A)f 3 - A mf -ai-5v,g]^+ivfVg= , (2) 

g®g = -n 2 l , (3) 

where A(R) is the vector potential of the static magnetic field, B = V x A, 



4 m /(p/,R;i) the mean-field order parameter matrix, and o"i(p/, R; e, t) is the 
impurity self-energy. The perturbation 5v(pf, R; t) includes the external electric 
field and the field of the charge fluctuations, 5p(R;t), induced by the external 
field. For convenience we describe the external electric field by a vector potential 
(5A(R; t) and the induced electric field by the electro-chemical potential <5</?(R; t). 
Hence in the Nambu-Keldysh matrix notation the perturbation has the form, 

5v = -%f5A(R;t)f 3 + e5ip(R;t)l , (4) 

and is assumed to be sufficiently small so that it can be treated in linear response 
theory. 

Equations (2) and (3) must be supplemented by self-consistency equations 
for the order parameter and the impurity self-energy. We use the weak-coupling 
gap equations, 

4£/( P/ ,R;t) = ^(nP/,P/)f(P/,R;M)) , (5) 

A« f (p f ,R;t) = , (6) 
and the impurity self-energy in Born approximation with isotropic scattering, 

ff,(R;M)^(j(p;,R;M)), (7) 

where f K is the off-diagonal part of the 2x2 Nambu matrix g K , and the Fermi 
surface average is defined by 

1 f rf2p / 

( - >= N~fJ • (8) 

4 The commutator [A,_B](g, is given by A Cg> B — B (g) A, were the noncommutative 
^-product is defined by A <g) B(e, t) = expi^ 8 ?"^ 9 ?) A(e, t)B(e, t) . 
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The materials parameters that enter the self-consistency equations are the pair- 
ing interaction, V(pf,p'j), the impurity scattering lifetime, r, in addition to the 

Fermi surface data p/ (Fermi surface), v/ (Fermi velocity), and Nf = J jjnp^vj]- 
In the linear response approximation one splits the propagator and the self- 
energies into an unperturbed part and a term of first order in the perturbation, 

g = go + Sg, A mf = A mf0 + SA mf , <7j = (T + 5&i , (9) 

and expands the transport equation and normalization condition through first 
order. In th order we obtain 

[(e + -v/ • A)f 3 - A mf0 - <to , g ] + iv f ■ Vg = , (10) 

go ® 50 = -7T 2 i , (11) 

and in 1 st order 

[(e+-v/-A)f 3 -4 m /o-d-o,^] (gl + iv/-V55 = [<54 m /+<5oi + <ft> , <?o]® , (12) 

go ® Sg + Sg®g = • (13) 

In order to close this system of equations one has to supplement the transport 
and normalization equations with the self-consistency equations of th and I s * 
order: 

^'/o(p/.R) = /7 £i(V<Pf>P'f)fo (P'/.B-Ie)), = o , (14) 

/+e c J, 
— (y(p / ,p' / )^(p^R;e,t)), 5^=0, (15) 

and 

<7 (R;e) = ^(j (p'/,R;f)) , (16) 

Sa i (R;e,t) = ^-(Sg(p' f ,B.-,e,t)) . (17) 

Finally, the electro-chemical potential, Sip, is determined by the condition 
of local charge neutrality (Gorkov & Kopnin 1975, Artemenko & Volkov 1979). 
This condition follows from the expansion of charge density to leading order in 
the quasiclassical expansion parameters. One obtains 5p(R;t) = 0, i.e. 

/j 
-^(TrSg K (p' f ,R;e,t)) = . (18) 
47T2 

The self-consistency equations (14)-(17) are of vital importance in the context 
of this paper. Equations (15) and (17) for the response of the quasiclassical self- 
energies are equivalent to vertex corrections in the Green's function response 
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theory. They guarantee that the quasiclassical theory does not violate funda- 
mental conservation laws. In particular, (16) and (17) imply charge conserva- 
tion in scattering processes, whereas (14) and (15) imply charge conservation 
in a particle-hole conversion process. Any charge which is lost or gained in a 
particle-hole conversion process is balanced by the corresponding gain or loss 
of condensate charge. It is the coupled quasiparticle dynamics and collective 
condensate dynamics which conserves charge in superconductors. Neglect of the 
dynamics of either component, or use of a non-conserving approximation for the 
coupling of quasiparticles and collective degrees of freedom leads to unphysical 
results. Condition (18) is a consequence of the long-range of the Coulomb re- 
pulsion. The Coulomb energy of a charged region of size £q and typical charge 
density eNfA is ~ e 2 N 2 A 2 t^, which should be compared with the condensation 
energy ~ NfA 2 ^. Thus, the cost in Coulomb energy is a factor (Ef/A) 2 larger 
than the condensation energy. This leads to a strong suppression of charge fluc- 
tuations, and the condition of local charge neutrality holds to very good accuracy 
for superconducting phenomena. 

Equations (10)- (18) constitute a complete set of equations for calculating the 
electromagnetic response of a vortex. The structure of a vortex in equilibrium 
is obtained from (10), (11), (14) and (16), and the linear response of the vortex 
to the perturbation 5A(R;t) follows from (12), (13), (15), (17) and (18). The 
currents induced by JA(R;i) can then be calculated directly from the Keldysh 
propagator Sg K via 



3 Structure and Spectrum of a Pancake Vortex in 
Equilibrium 

We consider an isolated pancake vortex in a strongly anisotropic, layered super- 
conductor, and model the Fermi surface of these systems by a cylinder of radius 
pj, and a Fermi velocity of constant magnitude, Vf, along the layers. We also 
assume isotropic (s-wave) pairing and isotropic impurity scattering. Thus, the 
materials parameters of the model superconductor are: T c , Vf, the 2D density of 
states Nf = pf/2-irvf, and the mean free path I — VfT. The model superconduc- 
tor is type II with a large Ginzburg-Landau parameter k 100, so the magnetic 
field is to good approximation constant in the region of the vortex core. 

We first present results for the equilibrium vortex. The order parameter, 
Z\ (R), is calculated by solving the transport equation (10), the gap equation 
(14), and the self-energy equation (16) self-consistently. These calculations are 
done at Matsubara energies (e — ► ie n = i(2n + l)nT). Details of the numerical 
schemes for solving the transport equation self-consistently are given elsewhere 
(Eschrig 1997). Charge conservation for the equilibrium vortex follows from the 
circular symmetry of the currents. Nevertheless, self-consistency of the equilib- 
rium vortex is important; the equilibrium self-energies (A , ctq) and propagators 




(19) 
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Fig. 1. Results of self-consistent calculations of the modulus of order parameter (a) 
and the current density (b) at T = 0.3T C for superconductors with electron mean free 
path £ — 10£o (clean) to i = |£o (dirty). R is the distance from the vortex center 
measured in units of the coherence length. 



(g ) are input quantities in the transport equation (12) for the linear response. 
Charge conservation in linear response is non-trivial and requires self-consistency 
of both the equilibrium solution and the solution in first order in the perturba- 
tion. 

Fig.l shows the order parameter and the current density in the vortex core 
of an equilibrium vortex for different impurity scattering rates r. As expected, 
scattering reduces the coherence length and thus the size of the core, and has a 
strong effect on the current density. Numerical results for the excitation spectrum 
of bound and continuum states at the vortex together with the corresponding 
spectral current densities are shown in Fig. 2. The local density of states (per 
spin) and the spectral current density are defined by 

N(K,e) = N f ^-Tr(r 3 (g R (p' f ,R,e)-g A (p' f ,R 7 e))) (20) 
j(R, e) = eNf^-Trfavftp'f) (g R (p' f ,K, e) g A (p' f , R, e)) ) . (21) 

The zero-energy bound state is remarkably broadened by impurity scattering. 
Its width is well approximated by the scattering rate l/(vf£), which is 0.63T C for 
£ = 10£o- The broadening of the bound states decreases with increasing energy. 
The results for the spectral current density (Fig. 2b) show that nearly all of the 
current density of the equilibrium vortex resides in the energy range of the bound 
states. This reflects the observation that the supercurrents in the vortex core are 
predominantly carried by the bound states (Bardeen et al. 1969; Rainer et al. 
1996). The physics of the core is dominated by the bound states, in particular 
also, as we will show, the response of the core to an electric field. 
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Fig. 2. Local density of states N(R, e) (a) and local spectral current density j x (R, e) (b) 
in the vortex core of a superconductor with £ — 10£o at T — 0.3T C . Results are shown 
for a series of spatial points on the y-axis, at distances 0, .257r£o, .57r£o, • • • , 47r£o from 
the vortex center. The thickest full line corresponds to the vortex center, and decreasing 
thickness indicates increasing distance from the center. Results for the outermost point 
(47r£o) are shown as dashed lines. The thin dotted lines show the density of states of 
the homogeneous superconductor (a) and the value of the bulk gap (b) respectively. 



4 Distribution Functions 



The set of formulas for calculating the quasiclassical linear response of a super- 
conductor is given in a compact notation in (9)-(18). In this section we transform 
these formulas into a more suitable form for analytical and numerical calcula- 
tions. The central differential equation of the quasiclassical response theory is 
the transport equation (12). It comprises 12 differential equations for the com- 
ponents of the three 2 x 2-Nambu matrices, Sg R ' A ' K . The number of differential 
equations can be reduced significantly by using general symmetry relations (see 
Rainer & Sauls, this volume) and the normalization conditions of the quasiclassi- 
cal theory. We focus on the simplifications which follow from the normalization 
equations (11) and (13). We use the projection operators introduced by She- 
lankov (1984), 

^-H'+i^-^-K 1 -^) ' (22) 
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Obviously, P RA + p RA — J. The algebra of the projection operators follows 
from the normalization conditions. 



pR,A pR,A 



pR,A pR,A 



\R,A 



0. 



(23) 



A key result is that the Nambu matrices Sg R ' A ' K can be expressed, with the 
help of Shclankov's projectors, in terms of 6 scalar distribution functions, Sj R ' A , 
5^ R ' A , Sx a and 8x a , each of which is a function of p/, R, e, t, and satisfies a 
scalar transport equation. The distribution functions are defined by 



and 



6g R ' A 
= ^2ni 



pR,A ( 



5i R > A 




R,A 



?R,A, 



( 




®P R - A 




o) 





(24) 



Sg a 

= -2iri 



P R < 



5x a 




P J 





0Sx a 



P A 



(25) 



where the anomalous response, 5g a , is defined in terms of 5g K , Sg R , Sg A by 

Sg K = Sg R <8> tanh(/3e/2) - tanh(/3e/2) <8> <55 A + <5.g a , (26) 

The transport equations for the various distribution functions follow from (10) 
and (12) and one finds (Eschrig 1997), 



ivj • V(?7 



R,A 



2eSj 

^R,A 



R.A 



+ (^ % ' /i A^ A - Z H > A ) ® <5 7 K < A + SY^ ® (4*^7^ + £^ A ) 
= -70^'^ ® <5^' A ® 70^'^+ ® 70^'^ 7o fl ' A ® ^ i? " 4 - (27) 



iv/ • V6^ R ' A - 2e^7 fl ' A 



+ (7 fl " 4 Z^" 4 - Z H > A ) ® 5 7 



#7 



7o 



(g) .W^" 4 - «5^ fl " 4 , (28) 



_^R,A ^ ^fl^A ^ , £ frfl.A ^ ~,fl,A ~R.A 



iv f ■ VSx a + id t Sx a + (^A R - S R ) ® &r a + fe a ® (Z\ A 7 A + S A ) 

= 7^ ® <5£ a ® 7^ - <5Z\ a ® 7^ - 7( f 5A a - 5S a , (29) 



iv f ■ V8x a - id t Sx a + (j R A R - S R ) ® 5x a + Sx a ® (A A -f A + S A ) 

= 7( f ® <5I7 a ® 7^ - <5^ a ® 7^ - 7( ? ® <5Z\ a - 5S a . (30) 
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We have used the following short-hand notation for the driving terms in the 
transport equations, which includes external potentials, perturbations and self- 
energies: 

e - ra ( £ B " A A R ' A \ 

--v/ • Af 3 + A mf0 + <t i0 ' = ( _^ R>A £ R . A j , (31) 

.r.a . / SS R ' A 8A R ' A \ _ fSU a SA a 



SA mf + 5ar + 5v={ _ sAR , A ^ R . A ) , 5af = { ^ ) . (32) 

The functions ^q' A and 7^ ' A in (27)- (30) are defined by the following convenient 
parameterization of the equilibrium propagators (Nagato et al. 1993, Schopohl 
et al. 1995): 

1 /-, R,A~R,A r, R,A \ 

~R,A _ . 1 / 1 - 7o 7o 2 7o \ (o^ 

After elimination of the time-dependence in (27)-(30) by Fourier transform one 
is left with four sets of ordinary differential equations along straight trajectories 
in i?-space. For given right-hand sides these equations are decoupled, and deter- 
mine the distribution functions 5~/ R ' A , S^f R ' A , 8x a and Sx a . On the other hand, 
the self-consistency conditions relate the right-hand sides of (27)- (30) to the so- 
lutions of (27)-(30). Hence, equations (27)-(30) may be considered either as a 
large system of linear differential equations of size six times the chosen number 
of trajectories, or as a self-consistency problem. We solved the self-consistency 
problem numerically using special algorithms for updating the right hand sides. 
Details of our numerical schemes for a self-consistent determination of the re- 
sponse functions are given elsewhere (Eschrig 1997). 



5 Response to an A.C. Electric Field 

We consider a pancake vortex in a layered s-wave superconductor and calculate 
the response of the electric current density, <5j(R;i), to a small homogeneous 
a.c. electric field, <5E w (t) = <5E exp(— iuot). The results presented in this section 
are obtained by solving numerically the quasiclassical transport equations (27)- 
(30) together with the self-consistency equations (15), (17), and the condition of 
local charge neutrality (18). The calculation gives the local conductivity tensor, 
<7jj(R, w), defined by 

5j i (R,t)=a ij (R,uj)6E?(t) . (34) 

Figures 3 and 4 show the time development of the current pattern induced by an 
oscillating electric field in x-direction with time dependence SE x (t) — SE cos(wt). 
Results are given for a medium range frequency (w = 1.5 A) and a low frequency 
(u> = 0.3 A). Two features should be emphasized. At medium and higher fre- 
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Fig. 3. Snapshots of the time-dependent current pattern in the core of a pancake 
vortex at successive times, t = 0, j^T, j^T , j^T , j^T, -^T (upper left pattern to lower 
right pattern) . The arrows show the current density induced by an a.c. electric field 
in x-direction, E x (t) = 8E cos(uii), of frequency to — 2n/T — 1.5A. The data are 
calculated in linear response approximation for a clean superconductor (I — 10£o) at 
T = 0.3T C . Current patterns in the second half-period, t = \T — j^T, are obtained 
from the patterns in the first half-period by reversing the directions of the currents. 
The distance between two neighboring points on the hexagonal grid is 0.257r£o- 



qucncies the current flow induced by the electric field is to good approximation 
uniform in space and phase shifted by tt/2 (non-dissipative currents). The phase 
shift of 7r/2 is the consequence of a predominantly imaginary conductivity at 
frequencies above A, as shown in Fig. 5. The current pattern at low frequencies 
(Fig. 4) is qualitatively different. At the vortex center the current is phase shifted 
by « 7r/4 in accordance with the conductivity at oj — 0.3Z\, which has about 
equal real (dissipative) and imaginary (non-dissipative) parts, whereas further 
away from the center the conductivity becomes more and more non-dissipative. 
Fig. 4 shows that the current flow at low frequencies is non-uniform. The dissi- 
pative currents exhibit a dipolar structure with enhanced currents at the vortex 
center, and back-flow currents away from the center. On the other hand, non- 
dissipative currents are approximately uniform. 

The frequency dependence of the local conductivity, a xx (R, u>) , for R along 
the x- and y-axis is shown in figures 5a, b. These figures include, for compar- 
ison, the Drude conductivity of the normal state and the conductivity of the 
homogeneous bulk superconductor. A significant feature of conductivity in the 
core is the strong increase of IZea at low frequencies. The conductivity is in 
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Fig. 4. The same as in fig. 3, but for a smaller external frequency, lo — 0.3A The 
length of the current vectors is scaled down by a factor 5 compared to those in fig. 3. 



this frequency range much larger than the normal state Drude conductivity and 
the exponentially small conductivity of the bulk s-wave superconductor. The 
real part of conductivity scales at low frequencies like 1/uj 2 . Its value at at the 
vortex center is 69.5e 2 Nfvj/A (outside the range of the figure) for w = 0.1 A. 
The enhancement of the dissipative part at low frequencies is a consequence of 
the coupling of quasiparticles and collective order parameter modes, and can- 
not be obtained from non-selfconsistent calculations. Fig. 5b shows that the real 
part of the conductivity becomes negative in the region of dipolar backfiow on 
the y-axis. This leads locally to a negative time averaged power absorption, 
< j • E > t < 0, and a corresponding gain in energy, which is compensated by 
the strongly enhanced dissipation of energy in the center of the vortex. The dis- 
sipative part of the local conductivity at a distance R from the vortex center 
exhibits pronounced maxima whose frequencies increase with increasing R and 
are given by 2 x the energy at the maxima in the local density of states shown 
in Fig. 2a. Hence, these features in the absorption spectrum must be identified 
as impurity assisted transitions between corresponding bound states at negative 
and positive energies. Impurities are required for breaking angular momentum 
conservation in these transitions. 

The applied electric field <$E(t) induces in the vortex core an internal field 
— V<fy>(i), which is of the same order as the applied field. Fig. 6 shows the total 
electric field, <5E tot (t) = £E w (t) — VS(p(t), in the vortex core. The induced field 
is at low frequencies of dipolar form, and oscillates out of phase (phase shift n/2) 
with the applied field. This dipolar field originates from small charge fluctuation 
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Fig. 5. Frequency dependence of the real part ('lower' curves) and imaginary part 
('upper' curves) of the local conductivity a xx for a superconductor with parameters 
T — 0.3T C , I — 10£o. For convenience, the conductivities are multiplied by uj. The 
full black curves give the conductivity at the vortex center, and the series of dashed 
lines with decreasing intensity the conductivity at increasing distance from the center. 
Fig. 5a presents data at points along the x-axis (in steps of (tt/4)£o), and Fig. 5b at 
points along the y-axis (in steps of (VStt/A)^)- The dashed grey lines show the normal 
state Drude conductivity, and the full grey lines the conductivity of the homogeneous 
superconductor. 



in the vortex core. At higher frequencies the dipolar field oscillates with a phase 
shift of « 7T, and screens part of the applied field. 

We finally discuss the role of self-consistency in our calculation. Our re- 
sults were obtained by iterating the self-consistency equations until the rel- 
ative error stabilized below < 10~ 10 . Fig. 7 compares the degree of violation 
of charge conservation in a non-selfconsistent calculation (no iteration) with 
the self-consistent result. We measure the degree of violation at position R by 
D(R) = maxi[(5p(R, t) + V • <5j(R, t)]. Charge conservation is obviously fulfilled 
if -D(R) = 0. The degree of violation at a point R is indicated in Fig.7a,b by 
the size of the filled circles around the grid points. The non-selfconsistent calcu- 
lation (Fig. 7a) results in a -D(R), which is much larger than the time derivative 
of the correct charge density, 5p(R, t) 7 obtained from a self-consistent calcu- 
lation with Sip(H, t) = instead of charge neutrality condition. In Fig. 7b we 
show the violation of charge conservation for our self-consistent calculation. The 
small remaining D(R) is here a consequence of the finite grid size used in our 
calculations. 
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Fig. 6. Local electrical field for a superconductor with parameters T = 0.3T C , i = 10£o, 
and an external frequency of u> = 0.3 A(T). Each field pattern is a snapshot at time 
t varying from to half of time periode in timesteps of periode. The external 
electrical field 5E w (t) = S~Ecos(cot) is maximal and points in positive ^-direction for 
the first picture (t = 0). In the first picture of the second line it is zero. The distance 
between two points in the grid corresponds to 0.257r£o- 

6 Discussion 

This paper presented self-consistent calculations of the electromagnetic response 
of a pancake vortex in a superconductor with finite but long mean free path. 
This complements previous calculations for perfectly clean systems, which were 
done self-consistently in the limit ui — > (Hsu 1993), and at finite frequencies 
without a self-consistent determination of the order parameter (Janko & Shore 
1992, Zhu et al. 1993). 

The frequency range of interest in our calculations is of the order of the 
gap frequency, hio = A. We have shown that at low frequencies (hui < 0.5A) 
the electromagnetic dissipation is strongly enhanced in the vortex cores above its 
normal state value, and that this effect is a consequence of the coupled dynamics 
of low-energy quasiparticlcs excitations bound to the vortex core and collective 
order parameter modes. The induced current density has at low frequencies a 
dipolc-like behaviour, which results from an oscillating motion of the vortex 
perpendicular to direction of the driving a. c. field. The response of the vortex in 
the intermediate frequency range, .5A ^ hui 5s 2A, is dominated by bound states 
in the vortex core. We find peaks in the local dissipation at twice the bound state 
energies. At higher frequencies, Huj > 2A, the conductivity approaches that of a 
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Fig. 7. Degree of violation of charge conservation by the dissipative current flow 
for a non-selfconsistent calculation (Fig. 7a), and a a self-consistent calculation 
(Fig.7b). The largest deviation in the non-selfconsistent calculation amounts to 
5p + V5j = 2.5e 2 NfVf8E" . The data are obtained for a superconductor with i — 10£o 
at u = 0AA(T) and T = 0.3T C . 

very clean homogeneous superconductor, which is in good approximation given 
by the non-dissipative response of an ideal conductor. 
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